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Abstract 

A phenomenological model was developed to explain quantitatively, with- 
out free parameters, the production of primary defects in silicon after particle 
irradiation, the kinetics of their evolution toward equilibrium and their in- 
fluence on detector parameters. The type of the projectile particle and its 
energy is considered in the evaluation of the concentration of primary defects. 
Vacancy-interstitial annihilation, interstitial migration to sinks, vacancy - im- 
purity complexes {VP, VO, V2O), and divacancy (V2) formation are taken 
into account in different irradiation conditions, for different concentrations of 
impurities in the semiconductor material, for 20 and °C. The model can 
be extended to include other vacancy and interstitial complexes. The density 
of the reverse current in the detector after irradiation is estimated. Com- 
parison with experimental measurements is performed. A special application 
considered in the paper is the modelled case of the behaviour of silicon de- 
tectors operating in the pion field estimated for the LHC accelerator, under 
continuum generation and annealing. 

PACS: 

29: Experimental methods and instrumentation for elementary-particle and 
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81: Materials science 
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1 Introduction 

The use of silicon detectors in high radiation environments, as to be expected in 
future high energy accelerators, poses severe problems due to changes in the prop- 
erties of the material, and consequently influences the performances of detectors. 

The incident particle, hadron or lepton, interacts with the electrons and with 
the nuclei of the semiconductor lattice. It losses its energy in several processes, 
which depend on the nature of the particle and on its energy. The effect of the 
interaction of the incident particle with the target atomic electrons is ionisation, 
and the characteristic quantity for this process is the energy loss or stopping power. 
The nuclear interaction between the incident particle and the lattice nuclei produces 
bulk defects and this phenomenon is studied in the present paper. As a result of 
this interaction, if the primary projectile is a particle, one or more light particles 
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arc formed, and usually one (or more) heavy recoil nuclei. This recoil nucleus 
has charge and mass numbers equal or lower than that of the medium. After 
this first interaction, the recoil nucleus or nuclei are displaced from the lattice 
positions into interstitials. Then, the primary knock-on nucleus, if its energy is large 
enough, can produce the displacement of a new nucleus, and the process continues 
as long as the energy of the colliding nucleus is higher than the threshold for atomic 
displacements. We denote these displacement defects, vacancies and interstitials, 
as primary defects, prior to any further rearrangement. In silicon these defects 
are essentially unstable and interact via migration, recombination, annihilation or 
produce other defects. 

As a consequence of the degradation to radiation of the semiconductor mate- 
rial, an increase of the reverse current due the reduction of the minority carrier 
lifetime, a reduction of the charge collection efficiency and a modification of the 
effective doping, due to the generation of trapping centres, are observed in the de- 
tector characteristics. In this paper, for the first time, a phenomenological model 
was developed to explain quantitatively, without free parameters, the mechanisms 
of production of the primary defects during particle irradiation, the kinetics of their 
evolution toward stable defects and equilibrium and the influence of the defects on 
detector parameters. The effects of the incident particle type, of its kinetic en- 
ergy and of the irradiation conditions on the concentration of defects are studied. 
Vacancy-interstitial annihilation, interstitial migration to sinks, vacancy - impu- 
rity complexes {VP, VO and V2O), and divacancy (V2) formation are considered 
in different irradiation conditions, for different concentrations of impurities in the 
semiconductor material and at different temperatures near room temperature. The 
model can be extended directly to include the effects of other mechanisms related to 
these, or other impurities in silicon, and of their interaction with the vacancy and/or 
interstitial. The density of the reverse current in the detector after irradiation is es- 
timated. Comparison with experimental published data of the time evolution of the 
concentration of defects is performed, as well as with measurements of the density 
of the leakage current. For different discrepancies, some explanations are suggested. 
A special application considered in the paper is the simulated case of the behaviour 
of silicon detectors operating in the pion field simulated for the future conditions 
at the new LHC (Large Hadron Collider) accelerator. 

2 Production of primary defects 

A point defect in a crystal is an entity that causes an interruption in the lattice 
periodicity. In this paper, the terminology and definitions in agreement with M. 
Lannoo and J. Bourgoin are used in relation to defects. 

The basic assumption of the present model is that vacancies and interstitials 
are produced in materials exposed to radiation in equal quantities, uniformly in the 
bulk of the sample. They are the primary radiation defects, being produced either 
by the incoming particle, or as a consequence of the subsequent collisions of the 
primary recoil in the lattice. 

The concentration of the primary radiation induced defects per unit fiuence 
(CPD) in the semiconductor material has been calculated using the explicit formula 
(see details, e.g. in references ||): 

CPD iE) = ^ fy(^) L{Ej,,)s^dn ^ ^^^^NIEL (E) (1) 

^ ^ 2Es^J ^\dn)^S^ 2Es^ ^ ' ^' 

where E is the kinetic energy of the incident particle, Nsi is the atomic density 
in silicon, is the silicon atomic number, Egi - the average threshold energy 
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for displacements in the semiconductor, Em - the recoil energy of the residual 
nucleus produced in interaction i, L{E]^i) - the Lindhard factor that describes the 
partition of the recoil energy between ionisation and displacements and {da /d£t)i 
- the differential cross section of the interaction between the incident particle and 
the nucleus of the lattice for the process or mechanism i, responsible in defect 
production. Na is Avogadro's number. The formula gives also the relation with 
the non ionising energy loss {NIEL). It is important to observe that there exists a 
proportionality between the CPD and NIEL only for monoelement materials. 

For CPD produced by pions, the pion - silicon interaction has been modelled 
01 and the energy dependencies of the Lindhard factors have been calculated in the 
frame of analytical approximations for different recoils in Si 

The concentration of primary defects produced by protons, neutrons, electrons 
and photons have been obtained from the NIEL. The calculations of Summers 
and co-workers for proton and electron NIEL in silicon from reference the 
calculations of proton, electron and photon NIEL of Van Ginneken Q as well as 
those of Ougouang for neutrons [|| have been considered. 

In Figure |l|. the dependence of the CPD on the particle kinetic energy is pre- 
sented: for pions, our calculations from reference |^ have been used; for protons, in 
the energy range 10^"^ 1 MeV the calculations of Summers, in the range 1 ~ 200 
MeV the average between those of Summers and of Van Ginneken, while in the 
range 200-^10000 MeV, Van Ginneken's. The curve for electrons uses, up to 1 MeV 
only the values from reference [|| , in the range 1 200 MeV an average between 
the values from references Q and , and after 200 MeV only from Q . The curves 
for photons and neutrons are calculated from Van Ginneken's and Ougouang's 
||] respectively. 

The main source of errors in the calculated concentration of defects comes from 
the modelling of the particle - nucleus interaction and from the number and quality 
of the experimental data available for these processes. Due to the important weight 
of annealing processes, as well as to their very short time scale, CPD is not a 
measurable physical quantity. 

In silicon, vacancies and interstitials are essentially unstable and interact via 
migration, recombination, annihilation or produce other defects. 

3 The kinetics of radiation induced defects 

In the frame of the model, equal concentrations of vacancies and interstitials are 
supposed to be produced by irradiation, in much greater concentrations than the 
corresponding thermal equilibrium values, characteristic to each temperature. Both 
the pre-existing defects and those produced by irradiation, as well as the impurities, 
are assumed to be randomly distributed in the solid. An important part of the 
vacancies and interstitials annihilate. The sample contains certain concentrations 
of impurities which can trap interstitials and vacancies respectively, and form stable 
defects. 

In the present paper, vacancy-interstitial annihilation, interstitial migration to 
sinks, divacancy and vacancy impurity complex formation (VP, VO, V2O), are 
considered. The mechanisms of formation of higher order defects involving vacancy 
and oxygen can be added, as well as the effects of other impurities, e.g. carbon. 

This picture could be described in terms of chemical reactions by the kinetic 
scheme: 

V + I ^annihilation (2) 
/ ^ sinks (3) 



3 




Figure 1: Energy dependence of the concentration of primary defects on unit fluence 
induced by protons, pions, electrons, photons and neutrons in silicon - see text for details. 

F + P g (4) 

VP is the E centre. 

V + O^^VO (5) 

VO is the A centre. 

V + V'^V2 (6) 

V + A^V20 (7) 

The bimolecular recombination law of interstitials and vacancies is supposed 
to be a valid approximation for the present discussion, because at the concentra- 
tions of vacancies of interest, only a small fraction of defects anneals by correlated 
annihilation if their distribution is random (see the discussion in reference |ll0| ). 

The multivacancy oxygen defects as, e.g. V3O, V2O2, V3O2, V3O3, are not 
considered in the model. 

The reaction constant A'l (corresponding to vacancy - interstitial annihilation) 
is determined by the diffusion coefficient of the interstitial atom to a substitutional 
trap: 

Ki = 30z/exp {-En/kBT) (8) 

where En is the activation energy of interstitial migration and v the vibrational 
frequency. The reaction constant in process (2) is proportional to the sink concen- 
tration a: 
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K2 = ai/A^cxp /fcsT) 



(9) 



with A the jump distance. 

Lee and Corbett argue that divacancies, vacancy-oxygen and divacancy- 
oxygen centres are equally probable below 350 °C; thus, K^, K^, Kj and Kg, that 
describe the formation of vacancy - impurity complexes and of divacancies, are 
determined by the activation energy of vacancy migration, Ei2, and are given by: 

Ks =K^ = Kt = Kg = 30iy exp (-^.a/fesT) (10) 

while /\4, Kq, Kg and Kiq are related to the activation energies of dissociation of 
the A, E, V2 and V2O centres respectively. 



K4 = 5i^cxp{-EA/kBT) 



Ke = 5vexp{~EE/kBT) 



(11) 
(12) 



Ks = 51/ exp {~Ev,)/kBT) (13) 

Xio = 5i.cxp {-Ev.o/ksT) (14) 

where Ea, Ee, Ev2 and EV2O are the dissociation energies of the A, E, V2 and V2O 
complexes respectively. 

G is the generation rate of vacancy-interstitial pairs, and is given by the product 
of CPD by the irradiation flux. Thermal generation is neglected, this approximation 
corresponding to high irradiation fluxes. 

In the simplifying hypothesis of random distribution of CPD for all particles, 
two different particles can produce the same generation rate for vacancy-interstitial 
pairs. 



G = [{GPD)^^^, ^ {E^)] ■ %art.a{El) = [(CPi?)p„rt.b {E2)] ■ %art.2iE2) (15) 

is fulfilled. 

Here, $ is the flux of particles (a) and (6) respectively, and Ei and E2 their 
corresponding kinetic energies. The system of coupled differential equations corre- 
sponding to the reaction scheme (2)- (7) cannot be solved analytically. 

The following values of the parameters have been used: En = 0.4 eV, Ei2 = 0.8 
eV, Ea = 1.4 eV, Eb = 1.1 eV, Ey^ = 1-3 eV, Ey^ = 1-6 eV, Ey^o = 10" Hz, A 
= 10^5 cm2, a = 10^° cm-2. 

Defect concentrations, as well as their time evolution, have been calculated solv- 
ing numerically the system of coupled differential equations. 

We would like to underline the specific importance of the irradiation and an- 
nealing history (initial material parameters, type of irradiation particles, energetic 
source spectra, flux, irradiation temperature, measurement temperature, tempera- 
ture and time between irradiation and measurement) on defect evolution. 

In Figures 2a d, the formation and time evolution of the divacancy, vacancy- 
oxygen, divacancy-oxygen, and vacancy-phosphorous is modelled in silicon contain- 
ing the initial concentrations of impurities: 10^'' P/cm'^ and 5x10^^ 0/cm'^, and 
irradiated with pions with about 200 MeV kinetic energy (corresponding to the in 
their maximum of CPD in the energetic distribution), at a total fluence of 10^^ 
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Figure 2: Time dependence of the concentrations of; a)V2, h)VO, c)V20 and A)VP, 
induced in silicon with 10^'' V /err? and 5x10^^ O/cm"^, irradiated with 200 MeV kinetic 
energy pions, at a total fluence of 10^"" pions/cm^ in different irradiation conditions - see 
text. 

pions/cm^, in different irradiation conditions. The instantaneous irradiation pro- 
cess is an ideal case where the total fluence is received by the material at time 
t = 0, and only the relaxation process is studied. This process supposes that the 
annealing effects are not present during irradiation. The second case considered is: 
the irradiation is performed in a single pulse for a time 2xf 0** seconds, followed by 
relaxation, and the third case is a continuum irradiation process with a generation 
rate of defects equal with 5x10^" pions/cm^/s. In these last two cases, annealing 
takes place during irradiation. During and after irradiation, the temperature is 
kept at 293K. As expected, after instantaneous irradiation the concentrations of 
defects are higher in respect with "gradual" irradiation. A special attention must 
be paid to the formation of VO and V2O centres, for which the same concentration 
is attained both in the case of continuous and instantaneous irradiation. 

The effect of the decrease of temperature, from 293 to 273 K during irradiation 
and annealing, is presented in Figures 3a and b. The material contains the same 
phosphorous and oxygen concentrations as in the modelled case presented in Figure 
2, and was irradiated with pions of 200 MeV kinetic energy, receiving continuously 
a fluence of 10^^ pions/cm^ in ten years, in accord to the pions simulated radiation 
field at LHC 

The increase of temperature increases the rate of all defect formation. In the 
case of VO and VP a plateau in the time dependencies is attained. Only for VP 
the plateau value is temperature sensitive. 

The rate of generation of primary defects in silicon influences the concentrations 
of all stable defects. In Figures 4a c, the time evolution of the concentrations 
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Figure 3: Time dependence of the concentrations of VO, VP, V2 and V2O induced in 
silicon with lO'^* P/cm"^ and 5x10^^^ O/cm'', irradiated with 200 MeV Itinetic energy pions 
a total fluence of 10^'' pions/cm'^ with the flux estimated for LHC, for 293 and 273 K. 



of VO, VP, V2 and V2O in silicon with lO" P/cm^ and 5x10^^ O/cm^, irradiated 
with pions of 200 IvIeV with the flux estimated for LHC, 10 and 100 times higher 
respectively, is presented. It can be observed that at the LHC generation rate for 
CPD, after 3x10'' seconds an equilibrium is established for the VP complex: its 
rate of formation equals its rate of dissociation - Fig. 4a. This time is shorter for 
higher generation rate, as can be observed in Figures 4b and 4c. The value of the 
plateau concentration for the vacancy- oxygen complex is attained after around the 
same time as the plateau for the concentration, in conditions of the LHC generation 
rate, and a shorter time, about 2x10^ sec. for a rate hundred times higher than the 
LHC one. For other defects, as divacancies and divacancy-oxygen, the processes to 
established the equilibrium are very slow. 

The formation of divacancy-oxygen is delayed in respect to vacancy oxygen, and 
for long exposure times, the same value for the concentration is obtained. 

The effect of oxygen in irradiated silicon has been a subject of intensive studies 
in remote past. In the last decade a lot of studies have been performed to investigate 
the influence of different impurities, especially oxygen and carbon, as possible ways 
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Figure 4: Time dependence of the concentrations of VO, VP, V2 and 1^20 induced in 
silicon with 10^"* P/cm^ and 510^^ 0/cm^, by continuous irradiation with 200 MeV kinetic 
energy pious with the flux: :a) estimated for LHC, b) 10 times the flux estimated for LHC, 
c) 100 times the flux estimated for LHC, at 293 K. 

to enhance the radiation hardness of sihcon for detectors in the future generation 
of experiments in high energy physics - see, e.g. references [ p^ |T^. These impu- 
rities added to the silicon bulk modify the formation of electrically active defects, 
thus controlling the macroscopic device parameters. If silicon is enriched in oxy- 
gen, the capture of radiation-generated vacancies is favoured by the production of 
the pseudo-acceptor complex vacancy-oxygen. Interstitial oxygen acts as a sink of 
vacancies, thus reducing the probability of formation of the divacancy related com- 
plexes, associated with deeper levels inside the gap. For this purpose, in the model, 
the effects of the initial oxygen concentration in silicon was studied. In Figures 5 
a, b, c, d the time dependencies of V2, VO, V2O and VP are presented, for silicon 
containing 10^^, 10^^, 10^^, and 10^^ atoms/cm^ initial oxygen concentrations. 

One can observe that vacancy-oxygen formation in oxygen enriched silicon is 
favoured in respect to the generation of V2, V2O and VP, confirming the considered 
hypothesis, so, for detector applications the leakage current is decreased. At high 
oxygen concentrations, the concentration of VO centres saturates starting from low 
fluences. 
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Figure 5: Effect of oxygen doping concentration on the time dependence of the concen- 
trations of: a)V2, h)VO, c)V20 and d)VP, induced in silicon with 10^* P/cm^ irradiated 
with 200 MeV kinetic energy pions at total fluence of 10^* pions/cm^ in one pulse. 



A difficulty in the comparison of model predictions with experimental data is the 
insufficient information in published papers regarding the characterisation of sili- 
con, and on the irradiation (flux, temperature during irradiation and measurement, 
irradiation time, time and temperature between irradiation and measurement) for 
most of the data. 

For electron irradiation, our simulations are in agreement with the measurements 
presented in reference p7[ | , where defect concentrations are presented as a function 
of the time after irradiation. In Figure 6, both measured and calculated depen- 
dencies of the VP and V2 concentrations are given. The irradiation was performed 
with 2.5 MeV electrons, up to a fluence of 10^^ cm~^. A good agreement can be 
observed for the concentration of VP^ see Figure 6a, while for the divacancy, Figure 
6b, the experimental data attain a plateau faster, and at smaller values than the 
calculations. The relative values are imposed by the arbitrary units of experimental 
data. 

A good agreement has also been obtained for hadron irradiation. For example, 
the sum of the calculated VP and V2 concentrations (8x10^^ cm^'^) induced in 
silicon by 5.67x10^'^ cm^^ 1 MeV neutrons, are in accord with the experimental 
value of 11.2x10^^ cm~^ reported in reference [|8|. 

4 Correlation with detector parameters 

It is well known that the dark current in a p — n junction is composed by three 
different terms: the diffusion current, caused by the diffusion of the minority charge 
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Figure 6: Time dependence of a,)VP and b)V2 concentrations after electron irradiation: 
points; experimental data from reference |ri] and dashed line: model calculations. 



carriers inside the depleted region; the generation current, created by the presence 
of lattice defects inside the bulk of the detector; and surface and perimetral currents, 
dependent on the environmental conditions of the surface and the perimeter of the 
diode. The appearance of the defects after irradiation corresponds therefore in an 
increase of the leakage current of the detector by its generational term. 

Inside the depleted zone, n,p << rii (n^ is the intrinsic free carrier concentra- 
tion), each defect with a bulk concentration iV^ causes a generation current per 
unit of volume of the form ll9ll: 



I = qU = q < vt > nt 



(TnCTpNT 



(16) 



where 7„ and gammap are degeneration factors, cr„ (cTp) are the cross sections for 
majority (minority) carriers of the trap, Ei ~ {Ec—Ev)/^ and < Wt > is the average 
between electron and hole thermal velocities. In the case of E and A centres and 
and V2 ) defects, the current concentration can be expressed in the simple 
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form: 



I = qU = q<Vt> nt^iVTet^'"^')/'^--^ (17) 

The primary effect in the recombination process is the change the charge state of 
the defect. The different charge states of the same deep centre may have different 
barriers for migration or for reacting with other centres. Thus, carrier capture 
can either enhance or retard defect migration or particular defect reactions. As a 
characteristics for detectors (as diode junction), the defect kinetics is dependent to 
the reverse - bias voltage during the irradiation [20[| . 

The comparison between theoretical and experimental generation current den- 
sities after irradiation shows a general accord between experiment and the model 
results for the lepton irradiation and large discrepancies for the hadron case. 

There could be several reasons for the observed discrepancies. 

The model hypothesis of defects distributed randomly in semiconductors exclude 
the possibility of cluster defects. For this case, other mechanisms of defect formation 
are necessary, which suppose different reaction rates and correlation between the 
constituent defects of the cluster. 

In the Shocklcy-Read-Hall model used for the calculation of the reverse current, 
each defect has one level in the gap, and the defect levels are uncoupled, thus the 
current is simply the sum of the contributions of different defects. In fact, the 
defects could have more levels, and charge states, as is the case of the divacancy, 
and also could be coupled, as in the case of clusters. As shown in the literature 
|2l| , [22I , both cases can produces modifications of the generation rate. 

Also the multivacancy oxygen defects as, e.g. V3O, V2O2, V3O2, V3O3, are not 
considered in the model. 

A model estimation of the time dependence of the leakage current, in conditions 
of continuous irradiations with pions of 200 MeV kinetic energy, in the conditions 
of the LHC ijl^, and at 293K is presented in Figure 7, for two concentrations of 
oxygen in silicon: 5x10^^ cm"'^ and 10^^ cm~^ respectively. As underlined before, 
oxygen incorporation in silicon has beneficial effects, decreasing the reverse current. 
This conclusion is valid in the hypothesis of random distribution of defects inside 
the depleted zone of the p-n junction. These values are probably underestimated. 



5 Summary 

A phenomenological model that describes silicon degradation due to irradiation, 
the kinetics of defects toward equilibrium, and the influence on the reverse current 
of detectors was developed. 

The production of primary defects (vacancies and interstitials) in the silicon bulk 
was considered in the frame of the Lindhard theory, and considering the peculiarities 
of the particle - silicon nuclei interaction. 

The mechanisms of formation of stable defects and their evolution toward equi- 
librium was modelled, and the concentrations of defects were calculated solving nu- 
merically the system of coupled differential equations for these processes. Vacancy- 
interstitial annihilation, interstitial migration to sinks, vacancy-impurities com- 
plexes {VP, VO and V2O), and divacancy formation were considered in different 
irradiation conditions, for different concentrations of impurities in the initial semi- 
conductor material and at different temperatures of irradiation. The calculated 
results suggest the importance of the conditions of irradiation, temperature and 
annealing history. The model supports the experimental studies performed to in- 
vestigate the influence of oxygen in the enhancement of the radiation hardness of 
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Figure 7: Time dependence of the reverse current after 200 MeV kinetic energy pions 
irradiation, with the rate estimated for LHC, at 293K, for silicon containing: a)5xl0^'' 
cm^'' and b)10^^ cm^'' oxygen. 

silicon for detectors. The VO defects in oxygen enriched silicon is favoured in re- 
spect to the other stable defects, so, for detector applications it is expected that the 
leakage current decreases after irradiation. The second result in the model is that 
at high oxygen concentrations, this defect saturates starting from low fluences. 

Most of the model calculations simulates some of the pion field estimated at 
the new LHC accelerator, where the silicon detector will operate under continuum 
generation and annealing. 

The density of the reverse current in detectors after irradiation is estimated, 
compared with experimental available data and for discrepancies some explanations 
are suggested. 
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